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ABSTRACT 

o 

Aims. We propose a method to be able to decide whether the planets of CoRoT-7 are moving on mutually inclined orbits in the order 
ofi> 10°. 

Methods. The extrasolar system CoRoT-7 is very special with respect to the closeness of the planets to the host star, which results 
in a fast dynamical development. It would therefore be possible to determine the inclination of the innermost planet CoRoT-7b with 
respect to the observer after an observation of at least three years from space with the satellite CoRoT with sufficient precision. 
' Different inclinations would cause different duration of the transit times of a planet in front of the star and would therefore give us a 

better knowledge of the architecture of this system. With the aid of numerical integrations and analytical estimations we checked how 
inclined orbits of additional planets would change the transit duration of CoRoT-7b 
i-C ■ Results. After 3 years of observations when an additional planet would be on a inclined orbit with respect to CoRoT-7b (/„„„„„/ > 10°) 

an increase of the order of minutes could be observed for the transit duration. 

O ' Key words, extrasolar planets, CoRoT-7, inclined orbits 
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1. Introduction 

(T) , Most of the more than 440 extrasolar planetary systems (=EPS£] are single planetary systems, but as of October 2010 there are 

■ some 50 multiplanetary systems, where between 2 and 7 planets are known to orb it their host stars. A very special case is the newly 
Q> | discovered system HD 10180 with seven planets orbiting an early G-type star (lLovis et al.ll2.010h . Certainly the small number of 

■ multiple EPS is a biased sample because it is highly improbable that there is just one planet around a star. From the theory of 
7 ~ I ' formation one expects that several planets may originate from the disk of gas and dust around a young star, like it is the case for the 
N J # | solar system. 

An interesting fact is that many of the EPS-planets seem to move on large eccentric orbits and consequently - when we expect 
CD more planets orbiting the star - strong perturbations will act on the planets and therefore the orbital elements may change signif- 
CD . icantly. Unfortunately our present methods to detect planets are just snapshots of the dynamical life time such that we don't have 
,— ' ' access to the dynamical evolution of a system. This is especially true for the detections via RV measuremenent where it is impossible 
| to determine the inclinations of the orbital planes with respect to the observer. Even via transit observations of one planet combined 
with the measurement of radial velocities it is impossible to determine mutual inclinations of the orbits of planets when only one 
- clearly the innermost - is transiting. The main problem to be solved would be to observe such an EPS for sufficiently long times 
to be able to see more than one snapshot in the dynamical evolution of the system. We are now in the same situation in astronomy 
as we had been during the last centuries when scientists wanted to determine the proper motion of stars or the orbital elements of 
double stars. Within the EPS most of the planets, especially the ones which we observe via transits, are relatively close to the host 
star and they have relatively large masses. So we might hope to see a signature of additional planets influencing the orbit of the 
inner planet which may result in a change of the transits times and/or the duration of the transit. Up to now it was impossible to 
make such precise measurements using ground-based observation but with the possibility of using satellites the situation changed. 
The CoRoT spacecraft (launched 2006) and the NASA Kepler mission (launched 2009) permit very accurate measurement of light 
curves of transiting planets in terms of photometric precision and timing of transits. The recent extension of the CoRoT mission for 
another three years is important with respect to the observation of Transit Time Variations (TTV) and also possible Transit Duration 
Variations (TDV) on a time line of several years. 
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2 The CoRoT space mission has been developed and is operated by CNES with the contribution of Austria, Belgium, Brazil. ESA, Germany 
and Spain 
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2. EPS with planets on inclined orbits 

The discovery of transiting planets with short periods makes searches for additional planets in an EPS via TTV and TDV feasible. 
There are a number of recent publications concerning the possible perturbations in such systems where outer planets with eccentric 
or even inclin ed orbits with respect to the innermost planet are present. 

Mardling (2010) was investigating three of such systems, namely HAT-P-13, HAT-P-7 and WASP-17. All three have orbital 
periods 2.2 days < P < 3.7 days and for HAT-P-13 an additional planets far away (P - 448 days) on a large eccentric orbit 
(<? = 0.666) is confirmed. The eccentricity of the inner planet, assuming coplanar orbits, may suffer from larger variations due to 
mutual inclinations of two planets, where the outer one is dominating the angular momentum. This effect of amplitude variations is 
increasing with increased mutual inclination. 

An important effect when two planets are moving on inclined orbits i s the Kozai resonance w hich has a stabilizing effect for 
orbits with 40° < i mu , < 60° having large eccentricities. In an article by Naga sawa et al.l (120081) the effect of coupling between 
mutual scattering, the Kozai mechanism and by tidal circularization was studied with the aid of orbital integrations. They developed 
a theory that a planet with small semimajor axis may not be formed by a type II migration because of planet-disc interaction, but 
through an excitation of the eccentricity because of the Kozai mechanism to values close to e — 1 such that the pericenter is close to 
the host star. Then the tidal friction may circularize the planet's orbit and as a consequence the authors predict that close-in-planets 
may have large inclinations up to retrograde orbits. 

iLibert & Tsiganisl (l2009al) treated the action of the Kozai resonance for five multiplanetary systems which are not in Mean 
Motion Resonances (=MMR), namely v Andromedae, HD 169830, HD 12661, HD 74156 and HD 155358. They found out by 
numerical studies in which they varied the unknown inclinations and also the nodal longitudes, that with the exception of HD 
155358 all of them could be in the Kozai resonance when the mutual inclination is larger than 45°. In another investigation by 
ILibert & Tsiganisl ([2009b) the role of mutually inclined orbits with respect to the stability of the orbits have been made in cases 
when the planets are in MMR. The authors have undertaken a parametric study varying masses and orbital parameters of the planets 
where they also took care of the migration rate and the rate of eccentriciy damping. It turned out that due to a ca pture process in the 
early phases of the system formation besides the 2/1 MMR (which was treated by |Th ommes & Lissauerll200l also in the 3/1, 4/1 
and 5/1 MMR the mutual inclinations may reach values as high as 70° when the eccentricity of one planet is larger than e > 0.4. Due 
to their results they say that our simulations show that inclination excitation is a common outcome, as long as eccentricity damping 
is not too strong. 

The case of TrEs-2b was expl ored by Scuderi et al.l (12 010); this is an EPS with a planet with a period P = 2.4 days. Using older 



tra nsit observations from 2006 bv lO'Donovan et a l. (2006) and comparing them to the ones of the Kepler mission recently published 
bv lGilliland et al.l (2010) they could not find significant changes in the orbital parameters since the discovery of that planet. 

I n another system, GJ 436, the transiting planet has a period of P = 2.6 days but a surprisingly large eccentricity of e = 
0.15. [Batygin etld] d2009h - using secular evolution of a two and also a three planet system describing GJ436 - succeeded in 
understanding this fact: the circularization has very long times scales when other planets on eccentric orbits are present which may 
even move in almost the same orbital plane as GJ436b. 

A different topic in astronomy, namely the observation of period variations in eclipsing binaries due to the presence of a third 
component, offers possi bilities of studyin g the so-called transit time variations (TTV) also in the analysis of transit photometry of 
extrasolar-planets (Bor kovits et alJuOlOl) . But also the duration of an eclipse was studied in the context of EB, and which is the 
consequence of an influence on the inclination of the two components with respect to the line of sights. 

The subject of this article, the extrasolar planetary system CoRoT-7 is hosting one planet with the very small period of P = 
20 hours and at least one other planet. We did numerical investigations for different sets of parameters which may cause a visible 
change in the transit observations especially with respect to the possible change of the inclination of the innermost planet CoRoT-7b. 



3. The extrasolar planetary system CoRoT-7 

The EPS CoRoT-7 contains on e transiting planet C oRoT-7b dLeger et al.ll2.Q09t) and one non transiting but confirmed planet CoRoT- 
7c detected by radial velocit y ( Oueloz et al. 2009). A reanalysis of RV data led to the assumption that there could be a third planet 
in this system dHatzes et alJl2010r) . CoRoT-7 is - up to now - a unique example for an EPS where the dynamical time scales are 
very short which is due to the closeness of the planets to the host star. The short periods of 0.85, 3.7 and 9 days for the three planets 
mean that the system is under very fast development: an integration time of a million years of CoRoT-7 in a numerical experiment 
is comparable to an integration of 100 million years for the solar system. Effects caused by mutual perturbations of the planets will 
be visible on short time scales. We list in Table Q] the orbital elements of the CoRoT-7 planets, showing the results obtained from 
two different analyses of the data. 



4. Investigation of mutually inclined planets in the CoRoT-7 system 

Former estimations showed that the regions where the CoRoT-7 planets are moving are very stable (e.g. lHatzes et aljfeOlOf) . We 
show that even inclinations of ij c b ~ 60° have no big influence on the eccentricity of the planets which would make the system 
unstable (see Fig. [3). The eccentricity changes are well inside the probable errors in the derived element eccentricity. The situation 
is completely different however for the inclination, because the duration of a transit is very sensitive to it. According to the published 
values of the inclination, respectively the errors (i = 80.1 ± 0.3° ), a shift of several tenth of degrees should be observable. In this 
section we will describe the numerical experiments and the results we obtained. 

It is particularly interesting whether one could determine how much the transiting planet would be influenced with respect to 
its inclination; this would have consequences for the transit time and the transit time duration. Therefore different test computations 
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Table 1. Orbital elements of the CoRoT-7 planets 



Parameters 



CoRoT-7b CoRoT-7c 
Data from 1 



CoRoT-7d 



mass [M Earth ] 
a[AU] 
P [days] 
e 

i [degree] 



4.8+0.79 
0.0172±0.00029 
0.853585±2.4xl0- 5 


80.1±0.3 

Data from ' 



8.39±0.89 
0.046 

3.698±0.003 




mass [M Earth ] 
a[AU] 
P [days] 



i [degree] 



6.9±1.4 

0.017 

0.853589 

±0.00059 

0.0 

80±0.3 3 



12.4±0.42 

0.045 

3.691 

±0.0036 

0.080±0.050 

<85 3 



16.7±0.42 

0.08 

9.021 

±0.019 

0.0±0.05 



1 http://exoplanet.eu 



Hatzes et al. (201C 
Oueloz et al.1 d2009h 



0.09 
0.08 
0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 



EXooCOOoooOOOoocOOoc 



xsCCOoooCXXXX 



10 



15 



20 



25 



30 




15 

time in 1000 years 



Fig. 1. Change in the perihelion and aphelion distances (upper graph) and inclinations (lower graph) for the CoRoT-7 planets with 
initial inclinations i-j c = ijj = 30° and i-ib =1° during 30,000 years (for colour figures see the online version of this journal). 



have been undertaken numerically (see Section l4~TI ) as well as using an analytical approach (see Section l4~2l for inclined planets in 
the CoRoT-7 system. 



4. 1 . Numerical results for the two planet case 

In this section we considered only the two confirmed planets in the CoRoT-7 system; an extension to the possible case of three 
planets can be found in Section l4~3l In Figs.[2][3]|4]and[5]we show the effect on the inclination of CoRoT-7b caused by a difference 
in the ascending node (Q72, = 20° and €l-i c = 140° = Run A, respectively vice versa = Run B). For each run, the following initial 
inclinations were used: in, — 1° and ij c = 30°, the time scale for the integration was 3.10 4 years. 

In Fig.|2]there is no visible difference between Run A and Run B over a time period of 30,000 years; note that CoRoT-7b reaches 
inclinations up to in, = 47°. A closer look at a shorter period of time (1,000 years) shows a slight time shift for the evolution of the 
inclination of CoRoT-7b (Fig.0. Zooming into the first 100 years of the integration it turns out that the time shift is caused by a 
different behaviour during the first ten years: while in, increases immediately in Run A, in Run B it decreases at first, but starts to 
increase after about ten years. From this point on the dynamical behaviour is qualitatively and quantitatively the same.: The thick 
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Fig. 2. Evolution of the inclination for CoRoT-7b and CoRoT-7c during 30,000 years of integration with initial inclinations iji, = l c 
and ij c = 30°. There is no visible difference between Runs A and B (for details see text). 




0.4 0.6 
time in 1000 years 

Fig. 3. Evolution of the inclination for CoRoT-7b and CoRoT-7c during the first 1,000 years of integration with initial inclinations 
i-jb = 1° and ij c = 30°. Run A and Run B are shown as thin and thick lines, respectively. 



and thin line in Fig. [4] show the same slope after the initial ten years. This time shift remains during the whole time which can be 
seen in Fig. [5] where we plotted the respective evolution of the inclinations after several thousand years. 
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Fig. 4. Detail of Fig. [3] Evolution of the inclination for CoRoT-7b and CoRoT-7c during the first 100 years of integration. Run A 
and Run B are shown as thin and thick lines, respectively. 
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Fig. 5. Evolution of the inclination for CoRoT-7b and CoRoT-7c for the last 400 years of the 30,000 year integration time. Run A 
and Run B are shown as solid and crossed lines, respectively. 
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Fig. 6. Secular change in the inclination of CoRoT-7b per year as a function of the relative position of the ascending node and the 
inclination of CoRoT-7c 



4.2. Analytical estimation 

We analytically estimate the change of the orbital elements in the two planet cas e (subscripts i and j)using the differential equations 
for the secular perturbation of the orbital elements <?, co, i, O from Stumpl (fl965t) : 

-j^ = -f [fi 8j ~ 8ifj) = di e j sin (w j - o>i) (4.1a) 

cos ii — = Cj sin z\ sin (O, - Q i) (4.1b) 
df v ; 

du>i e\ i \ 

—Ci-dj — cos [ojj — oji) (4.1c) 

df ei v ; 

d£2, ( sin ij 



df ' I sin;. 



1 cosfOi-O.H (4.1d) 

sin ti V J / 



In the case of a heavy central body (the star) and two planets, the right side of equations 14.11 can be assumed to be constant 
ignoring the mutual periodic perturbations. The orbital elements <?, o>, i, Q. can then be described as slowly increasing or decreasing 
functions of time. The differential equation for the inclination ^, lbl can thus be written as 

/ \ sin/,- 

Mi = -a sin (Q, - Q/) - (4.2) 

V ' COS/; 

depending only on the current values of /,, ij, O, and Clj. The c, are described by 

c i =^mjn i bf 2) (4.3) 

where is the Laplace coefficient of the first order. The b\ ^ depend on the ratio of the orbital distance of the two planets a 
and are calculated as follows: 

„ = 2 _ r cos^y 

* "Jo (1 -20-cosy + a 2 )- 5 

with a = 2i. 

a, 

Thus, for a given set of /,, ij, Q, and it is possible to estimate the change Ai using 14.21 Figs.|6]and|2]show the secular change 
in inclination of CoRoT-7b in degrees per year as a function of the relative position of the ascending node AQ = £lj - £2, and the 
inclination of the other planet. If we take into account only the perturbation of planet CoRoT-7c on CoRoT-7b, the change in its 
inclination occurs at a rate of ±0.2° per year at a separation of the ascending nodes of +90° . Considering only the perturbation of 
CoRoT-7d, the change is much smaller, namely ±0.06° per year (see Figs.|6]and|7]). The initial value for the inclination of CoRoT-7b 
was assumed to be 1° in both calculations. 

The results from this analytical approach fit the data from numerical integration very well, as can be seen in the following 
graphs: Figs.[8](analytical model) and[9](numerical model) show the change in A/7,1, after 3 years as a function of Q.j c and Q.jd, again 
with an initial inclination in, — 1° . 



4.3. Extension to the three planet case 

Since the data analysis of lHatzes et al.l d2010t) lead to the conclusion that (at least) three planets exist in the extrasolar planetary sys- 
tem CoRoT-7, we investigated the combined influence of different initial conditions for CoRoT-7c and CoRoT-7d on the inclination 
of the transiting planet. Here we ignored the possible differences in the ascending nodes of the three planets and set all to the same 
value. This assumption may be somewhat artificial but the parameter space is so large that this is the only way to estimate possible 
changes of the ij c . We varied the inclinations of both outer planets from normal prograde to retrograde orbits (0° < / < 180°). 
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Fig. 8. Change in the inclination of CoRoT-7b after 3 years as a function of the relative position of the ascending node of CoRoT-7b 
and CoRoT-7c using equation l4.2l 



- After 1 year: The results of the integration for such a short time interval (Fig.fTOll show that only for large inclinations of CoRoT- 
7c (i ~ 60°) CoRoT-7b would suffer from changes in the order of 1° without a big influence of the inclination of CoRoT-7d. On 
the other hand the outermost planet would have a significant influence only if its orbital plane was almost perpendicular to the 
plane of CoRoT-7b. This is also clear from the results of analytical computations (see Figs.|6]and|7]). 

- After 3 years: During the lifetime of the CoRoT satellite the change in the inclination of CoRoT-7b would be observable within 
the given error bars, when both outer planets would have inclinations in the order of several tenths of degrees (Fig.fTTTi. 

- After 1 years: It is evident from Fig.[T2]that for such a long time interval of observations the changes could be significant even 
for only slightly inclined orbits of one or both of the outer planets. 

The cross-shaped white area (for retrograde orbits) of the parameter space would lead to unstable orbits even in short time scales. 



5. Transit duration 



Sever al formulas have been published to calculate the duration of transits (e.g. I Kipping! l2Q08b ISeager & Mallen-Ornelasl 
2003b. To determine the transit duration of CoRoT-7b for different apparent inclinations we applied the following formula of 
Seager & Mallen-OrnelaJ ([2003J) 



arcsin 



A' 

a 



( 1 + f) 2 -(f COBi ?) 



l/2\ 



(5.1) 



being valid only for circular orbits. Note that in our calculations we always consider the total transit duration, i.e. including ingress 
and egress of the planetary disk. 

The numerical integrations of the CoRoT-7 system for ten years (see Fig.[T2l resulted in a change in the inclination of CoRoT-7b 
of up to 3° . Table|2]shows the corresponding change in the transit duration. 

The total duration of the transit, as calculated with Formula [5.1l - which in this case exactly des cribes the transit geometry - is 
66.4 minutes. This is in very good agreement with the transit duration published bv lLeger et al.l d2009l) of 1.125 + 0.05h = 67.5 + 3 
min. If the apparent inclination of CoRoT-7b decreases by a value of 3° , the transit duration decreases by approx. 32 minutes, 
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Fig. 9. Change in the inclination of CoRoT-7b after 3 years as a function of the relative position of the ascending node of CoRoT-7b 
and CoRoT-7c using results from numerical integrations 



CoRoT-7b: inclination after 1 years (d-con:0.5) 




inclination of CoRoT-7c 



Fig. 10. Change in the inclination of CoRoT-7b per year as a function of the relative position of the ascending node of CoRoT-7b 
and CoRoT-7c using results from numerical integrations. 



whereas if the inclination increases by 3° the transit duration is approx. 16 minutes longer (see Table |2j. This could be easily 
observed by the CoRoT satellite. Our numerical calculations suggest - depending on the exact conditions - that a 1° change in the 
inclination of CoRoT-7b could be expected after three years, and a 3° change should be possible over the course of 10 years (see 
Section l4~3l l. One a side note, if the apparent inclination of CoRoT-7b is less than 76.2° the transit cannot be observed any more. 

Since CoRoT-7 is a multiplanetary system with at least two planets - as already pointed out dOueloz et al.ll2009t lHatzes et alJ 
2010; see Table [D - we also investigated the possibility of CoRoT-7c or CoRoT-7d becoming a transiting planet. Taking the or- 
bital plane of CoRoT-7b to be the reference plane of the system, we calculated the transit possibility as a function of inclination 
and longitude of the ascending node in reference to CoRoT-7b. These calculations could not be done using formula I5T1 which 
takes into account only the apparent inclination, but not the longitude of the ascending node. This is an important orbital param- 
eter for determining the transit probability. We adopted a semi-analytical approach using the computer algebra software Wolfram 
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CoRoT-7b: inclination after 3 years (d-con:0.5) 




inclination of CoRoT-7c 



Fig. 11. Change in the inclination of CoRoT-7b per year as a function of the relative position of the ascending node of CoRoT-7b 
and CoRoT-7c using results from numerical integrations. 



CoRoT-7b: inclination after 10 years (d-con:2) 




inclination of CoRoT-7c 



Fig. 12. Change in the inclination of CoRoT-7b per year as a function of the relative position of the ascending node of CoRoT-7b 
and CoRoT-7c using results from numerical integrations. 



Mathematica™ . The planetary orbit is represented by a 3D parametric equation. The full representation of the orbit is obtained by 
combining a parametrization of the position of the orbital plane according to the Keplerian elements e, a, i, u> and Q, with a rotation 
of the system according to the point of view of the observer. This approach allows a consistent choice of the reference plane for 
multiplanetary systems, and accounts for the apparent inclination the system is observed under. We use an extended projection of 
the stellar disk in direction of the observer, represented by a cylindrical parametrization with a diameter of R — R* + Rp, so that we 
get the total transit duration. The intersection of the orbit with the projection of the disk yields two positions t\ and ?2 on the orbital 
curve, which correspond to the starting and ending points of the transit. For a circular orbit, which we assume for the planets of the 
CoRoT-7 system, the parameter t of the orbital curve is equivalent to the mean motion of the planet. In order to obtain the amount of 
time passing between the planet entering and exiting the stellar disk we calculate the length I of the arc of the orbital curve between 
the two points t \ and ti 
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i[degree] 11 A 11.6 78.1 78.6 79.1 79.6 80.1 1 80.6 81.1 81.6 82.1 82.6 83.1 

At, mns [min] 34.5 42.3 48.6 53.9 58.6 62.7 66.4 69.7 72.7 75.4 77.8 80.1 82.1 

1 currently adopted value of inclination (http://exoplanet.eul 
Table 2. Change in the transit duration of CoRoT-7b. 




-SO -<50 - 40 - 20 20 40 GO SO 
i D 

Fig. 13. Duration of a possible transit of CoRoT-7c as a function of inclination and longitude of the node. The transit duration is 
given in minutes and indicated by color coding and contour lines. 




Fig. 14. Detail of Fig. [T3l with labeled contour lines. 



Z = 2na- - (5.2) 

360° V ' 

The ratio of the orbital length during the transit to the circumference of the orbit L corresponds to the ratio of the duration of the 
transit to the planetary period P. The transit duration can then be written as 

t Transit = y P- (5.3) 

For investigating the possibility of CoRoT-7c or CoRoT-7d becoming a transiting planet, we had to adopt a value for the 
planetary radii which are currently not known. CoRoT-7c has a minimum mass of 8.4 M Earth, positioning it right between Earth and 
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Neptune. So we assumed the density of CoRoT-7c to be the mean value of the densities of Earth and Neptune, yielding a radius 
of 0.22 Rjupiter- CoRoT-7d has a minimum mass very similar to the mass of Neptune (16.5M£ nrt /,), so it was assumed to have the 
same density, corresponding to a radius of Q3ARj up i, er . Figures [13] and [14] show the duration of a transit of CoRoT-7c as a function 
of its inclination and the position of the ascending node. The reference plane is fixed to the plane of CoRoT-7b and the ascending 
node of CoRoT-7b is assumed to be in the direction of the observer. Thus, the value of the inclination as shown in Fig.[13]is not the 
apparent inclination as it would appear to an observer, but the in-system inclination with respect to the plane of CoRoT-7b. Fig.[T4lis 
showing a enlarged detail of Fig. [13] with labeled contours of equal transit duration. The same calculation was done for CoRoT-7d, 
and qualitatively the same behaviour was found. Due to the similarity of the plots - the only difference being in the scale - we 
refrain from showing this figure. The maximum transit duration is 153 minutes for CoRoT-7c, as can be seen in Figs. [T31 and [141 
and 217 minutes for CoRoT-7d. 

6. Conclusions 

In this investigation we estimated the TDV caused by inclined additional planets which reaches - although relatively small - values 
of some degrees and is therefore detectable with the CoRoT satellite. This possibility stems from the fact that the CoRoT-7 system 
is a rapidly evolving extrasolar system with very close-in planets. 

We did our study using numerical integrations for the long term development of the system, as well as an analytical approach 
for short time intervals in the order of years. It turned out the system is quite stable even in the 3 planet model. A quantitatively new 
result is the dependence on the difference in the ascending node on the short time development of CoRoT-7b. The small phase shift 
in the dynamical development is not important for the qualitative behaviour for long time, but essential for the short time evolution 
of the orbit of CoRoT-7b. 

In our determination of the duration of the transit caused by the change of the inclination of CoRoT-7b we also discussed what 
kind of orbits of the outer planet(s) would lead to a transit of these planets. 

Despite the constraints given by the incompleteness of the data derived from observations the main conclusion of our study is 
that after three years of observation of the EPS CoRoT-7 we would be able to determine - via transit observations from space of 
the CoRoT satellite and additional ground based RV - whether the two (three planets) are on mutually inclined orbits or whether 
they have just small inclinations of the order o fj < 10°. More work has to be done for a more detailed analysis with numerical 
studies but also with analytical approaches (e.g. Borkovits et al. 2003, 2007) comparable to those which have been undertaken for 
the investigation of the change in occultations of eclipsing binaries caused by an additional star. 
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